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Abstract 

It is shown that the large angle MSW solution of the solar neutrino problem with 
a bi-maximal neutrino mixing matrix implies an energy-independent suppression 
of the solar u e flux. The present solar neutrino data exclude this solution of the 
solar neutrino problem at 99.6% CL. 
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The possibility that the neutrino mixing matrix U has the bi-maximal mixing form 
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has attracted a large attention [Q] after the presentation at the Neutrino '98 conference of 
the Super-Kamiokande evidence in favor of atmospheric neutrino oscillations with large 
mixing 0. 

Neutrino bi-maximal mixing is capable of explaining in a elegant way the atmospheric 
neutrino anomaly 0, ||, £|, ||, through — > v T oscillations due to[| Am^ ~ 10~ 3 eV 2 and 
the solar neutrino problem (SNP) || 0, ||, || [K| through z/ e — > z/^, i/ T oscillations in 
vacuum due to Am^ ~ Kr 10 eV 2 [|TTJ. 

As noted in [pT2| , the results of the recent analysis of solar neutrino data presented 
13 ] seem to implyQ that neutrino bi-maximal mixing may be also compatible at 



111 



CL with the large mixing angle (LMA) MSW [|T4f solution of the SNP |]| (see Fig. 2 of 
0)- 



Here I would like to notice that this conclusion seems to be in contradiction with the 
exclusion at 99.8% CL of an energy-independent suppression of the solar v e flux presented 
in the same paper |Tj| (see section IV. D). 



The reason of this incompatibility is that bi-maximal mixing with the Am^ ~ 10~ 5 — 
10~ 4 eV 2 corresponding to the LMA solution of the SNP implies an energy-independent 
suppression by a factor 1/2 of the solar v e flux. 



This can be seen following the simple reasoning presented in |IE |. The mixing of the 
neutrino states in vacuum is given by (see, for example, [[nfl) 

W) = U akWk) (a = e,/i,r), (2) 

fc=l,2,3 

where the states |z/ a ) (a = e,fi,r) describe neutrinos produced in weak interaction pro- 
cesses and the states \uk) {k = 1,2,3) describe neutrinos with definite masses m^. 

In the bi-maximal mixing scenario the numbering of the massive neutrinos is the 
usual one, i.e. such that mi < vri2 < ms, and Am^ ~ lCT 3 eV for the solution of the 
atmospheric neutrino anomaly. If Am^ ~ 1CT 5 — 1CT 4 eV 2 for the LMA solution of the 
SNP, we have Amj 2 ~ Am^ ~ 10 _3 eV 2 . 

Solar neutrinos have energy E ~ 1 MeV and the ratio Am^/£ ~ Am| 2 /E ~ 10~ 9 eV 
is much larger than the matter induced potential V < 10~ n eV in the interior of the 
sun. Hence, the evolution equation of the heaviest massive neutrino z/ 3 is decoupled from 
that of the two light neutrinos v x and v 2 (see, for example, |]TS|). Taking also in account 
that in the case of bi-maximal mixing f/ e3 = 0, one can see that an electron neutrino 



1 Am^- = m^, — rrij is the difference between the squared masses of the two massive neutrinos and 
•j. In the bi-maximal mixing scenario there are three massive neutrinos, v\, and ^3. 

2 I want to emphasize from the beginning that 1 do not want to criticize at all the beautiful paper 
L3|. I am only concerned with the interpretation of its results. 
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is created in the core of the sun as a superposition of the two light mass eigenstates v\ 
and v 2 and, whatever happens during his propagation in the interior of the sun, its state 
when it emerges from the surface of the sun is a linear combination of \v\) and \i>2): 



W)s = ^2 a k Wk) , 



(3) 



fc=l,2 



with 
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Since the massive neutrino states \v k ) propagate as plane waves, the state describing the 
neutrino detected on the Earth is 



\v)e = ^ Qfce 

fc=l,2 



-iE k L 



Vk) 



(5) 



where L is the distance from the surface of the Sun to the detector on the Earth. The 
survival probability of solar electron neutrinos is then given by P^ 



sun 



\{VeV>)i 



a k e~ lEkL (is e \v k ) 



fc=l,2 



a k e~ lEkL U ek 



fe=l,2 



(6) 



Taking now into account the explicit values U e \ = l/v^2 and U e 2 = —l/y/2 in the case of 
bi-maximal mixing and the fact that the neutrinos are extremely relativistic, we have 



a\ — a 2 exp 



IE 



(7) 



In the case of the LMA solution of the SNP Am^ ~ 



10~ 5 — 10~ 4 eV and the survival 
probability (0) oscillates with an oscillation length AwE / Am^ ~ 10 7 cm that is about 
one million times smaller than the Sun-Earth distance. Hence, the oscillations are not 
observable on the Earth because of averaging over the energy spectrum and only the 
average probability 
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• > = §(N 2 + |02| 2 ) 



is observable. We have obtained the announced result: the LMA solution of the SNP in 
the bi-maximal mixing scenario implies an energy-independent suppression of the solar 
u e flux of a factor 1/2. 

Therefore, we have the apparent paradox that an energy- independent suppression of 
the solar v e flux seems to be allowed at 99% CL by Fig. 2 of Ref. |13|] and is excluded 
at 99.8% CL in Section IV.D of the same paper. Notice that the two conclusions are 
based on the same set of data and the same theoretical calculation of the neutrino flux 
produced by thermonuclear reactions in the core of the sun 

The fact that the two cases refer to the same physical situation, i.e. an energy- 
independent suppression of the solar v e flux, is also shown by the \ 2 calculated in the 
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two cases. The \ 2 °f the right border of the LMA region)] in Fig. 2 of Ref. is 
4.3 + 9.2 = 13.5, whereas the \ 2 calculated in Section IV.D of the same paper for an 
energy-independent suppression of the solar v e flux by a factor 0.48 is 12.0. Since this is 
the best fit for an energy-independent suppression of the solar v e flux, a value of \ 2 = 13.5 
for a suppression factor 0.5 looks plausible. 

The solution of the apparent paradox explained above lies in a correct statistical 
interpretation of the allowed LMA region in Fig. 2 of Ref. / [i^/ and of the exclusion in 
Section IV.D of the same paper. The two cases have different statistical meanings. 



The allowed regions in Fig. 2 of Ref. [y| are obtained under the assumption that the 



neutrino masses and mixing parameters are not known. In this case a general neutrino 
oscillation formula is used in the fit, with the neutrino masses and mixing angles consid- 
ered as free parameters. The best fit in the LMA region happens to have a Xmin = 4-3, 
which corresponds to a CL of 3.8% with 1 DOF. Hence, a LMA solution is allowed at 
3.8% CL and one can draw a 99% CL region corresponding to the parameters that have 
X 2 < Xln + 9-2. 

The statistical analysis discussed in Section IV.D of Ref. [O] assumes that the solar 



v e flux is suppressed by a constant factor that is the free parameter to be determined by 
the fit. It happens that the best fit has Xmm = 12.0, which corresponds to a CL of 0.2% 
with 2 DOF. Hence, the hypothesis is excluded at 99.8% CL and no allowed region of the 
free parameter can be drawn. 

Since the two statistical analyses start from different assumptions, it is clear that 
they answer different questions and their conclusions cannot be compared. Moreover, 
it is important to notice that the test of the maximal mixing scenario with Am^ ~ 
10 -5 — 10~ 4 eV 2 does not correspond to either of the two statistical analyses. Indeed, 
if this scenario is assumed, we know that the solar v e flux is suppressed by an energy- 
independent factor 0.5 and there is no parameter to fit. Hence we test the hypothesis 
under consideration on the basis of its \ 2 . The \ 2 ~ 13.5 indicated by Fig. 2 of Ref. |i3|| 
implies a CL of 0.4% with 3 DOF. Therefore, the hypothesis is rejected at 99.6% CL. 

Notice that this exclusion is based only on the values of the elements U e i, U e 2 and 
[/ e 3 of the neutrino mixing matrix. This means that also other types of neutrino mixing 



matrix, as those discussed in ||20|| , are incompatible with the LMA solution of the SNP. 

In conclusion, I would like to emphasize that the allowed regions of the neutrino 
oscillation parameters calculated in the usual way (i.e. as Fig. 2 of Ref. |13[) cannot be 



used to test a definite model (as the bi-maximal mixing model) because they have been 
obtained under different assumptions^. In order to obtain allowed regions appropriate for 
model testing one must use the procedure described in |21j, E?2| , i-e. one must consider 



each point of the parameter space as a model and perform a goodness of fit testing with 
it. I think that it would be very useful if both types of allowed regions will be presented 
in future papers. 



! If U e3 = 0, we have sin 2 2d = -i\U el \ 2 \U, 



< 2 



2 



(see and sin 2$ = 1 corresponds to \U e i \ = \U e 2\ = 

l/\/2, as in the bi-maximal mixing matrix (|l). 

4 They are useful if one wants to know the allowed range of the mixing parameters for other purposes. 
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